Standing shocks in the inner slow solar win 
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We examine whether the flow tube along the edge of a coronal streamer supports standing 
shocks in the inner slow wind by solving an isothermal wind model in terms of the Lambert W 
function. We show that solutions with standing shocks do exist, and they exist in a broad area in 
the parameter space characterizing the wind temperature and flow tube. In particular, streamers 
with cusps located at a heliocentric distance > 3.2i?0 can readily support discontinuous slow winds 
with temperatures barely higher than 1 MK. 
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It was proven possible that the quasi-steady solar wind 
may not be continuous but involve standing shocks in the 
near-Sun region First pointed out 30 years ago[l|, 
the existence of standing shocks depends critically on the 
existence of multiple critical points (CPs). These can 
arise due to either momentum addition or rapid tube 
expansion near the base. Time-dependent simulations 
showed that whether the system adopts a continuous or 
a discontinuous solution depends on the detailed manner 
the tube geometry is varied0,Q, or how the momentum 
addition is appliedjl, [1]. Existing studies on standing 
shocks were exclusively on the flow rooted in the inte- 
rior of coronal holes. However, little is known about 
whether the flow tubes bordering bright streamer hel- 
mets can support standing shocks as well. This region 
is important, however, since it is where the slow wind 
likely originates Q. Here the tube expansion is distinct 
from the coronal-hole one, with the tube likely to expe- 
rience a dramatic expansion around the streamer cusp 
(see Fig. 4, the current-sheet case in@). This letter is 
intended to answer: Are standing shocks allowed by this 
geometry? 

To isolate the geometrical effect, we will use a simple 
isothermal model. Let T and Vr denote the solar wind 
temperature and radial speed, respectively. The isother- 
mal sound speed is then Cg — \/2k^Tjrn^, where ks is 
the Boltzmann constant, and nip the proton mass. The 
Mach number M = Vr/cs is governed bvll2| 
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where y = t/Rq, with Rq the solar radius and r the 
heliocentric distance. Moreover, a = u/Rq is the non- 
dimensionalized tube cross-section a. And o is related to 
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the expansion factor / by a(r) — f{r)r'^. Furthermore, 
A = gQRQ/cj. where gQ is the surface gravitational accel- 
eration. Evidently A measures the relative importance 
of the gravitational force and pressure gradient force. 
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FIG. 1: Expansion factor / for the streamer geometry vs. 
heliocentric distance r. Please see text for the meaning of 
.foo, /a/, rc, and 5, and what the diamonds refer to. 



The streamer geometry is parameterized as 
1 



G{r-rcJ)-G{Re;rc,S) 
l-G{RQ;rc,5) ' 



ihl - 1) 

r < rc, 
foo + (/m - foo) G{r; rc,S), 
r > rc, 



(2) 



where G{x;xo,6) = exp 



is a Gaussian. 



{x-xo) /S^ 

Figure [T] illustrates the r-distribution of /. Obviously 
foo represents the value at large distances, and /m is 
the maximum attained at rc, the heliocentric distance 
of the streamer cusp. Moreover, S describes how rapid 
f]\i is approached. For rc, we adopt values between 2.4 
and 3.6i?0, compatible with LASCO C2 images. The 
ranges for foo, /m and 6 are [2, 10], [6, 22], and [0.4, 1]Rq, 
respectively. As direct measurements of the coronal mag- 
netic field remain largely unavailable, some model field 
is used to guide our choice. The / profile with the base 
values (/oo — 6, /m — 14, and 6 — 0.7Rq) is close to the 
diamonds in Figll] which correspond to / along the tube 
at the streamer edge in a current sheet model, given in 
Fig.4b ofd (the one labeled 27°). 
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Given the temperature T and an /(r), the right hand 
side (RHS) of Eq.([T]) can be readily evaluated and de- 
termines whether solutions with standing shocks are al- 
lowed. To explain this, we note that any root of RHS = 
corresponds to a critical point (CP), which is cither a lo- 
cal extreme (dM/dy = 0, M ^ 1) or a sonic point (SP) 
{dM/dy y^O,M = 1, denoted by the subscript S). Shock 
solutions are known to exist only when there are multiple 
CPs, and were usually constructed by carefully examin- 
ing the solution topology. Here we present a new method 
based on a recent study which shows that a transonic so- 
lution to Eq. (HI is expressible in terms of the Lambert W 
function W{xy 10] 



T=1.2 MK 



T=1.3 MK 



-Woi-Diy)), l<y<ys, 
-W-,{-D{y)), y>ys, 
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Only two things about W{x) need to be known in the 
present context: first, a real-valued W{x) can be defined 
only for a; > — 1/e (note that D is positive definite); sec- 
ond, W{x) has two branches for —1/e < a; < 0, and they 
obey — 1 < Wo < and W-i < — 1. The mathematical 
details can be found in [11]. In practice, we evaluate Wq 
and W^i via Eq.(5.9) there. 

If only one CP exists, it is naturally the SP, and Eq.Q 
describes the only possible transonic solution for which 
M increases monotonically with r. This is the case con- 
sidered in[l3|, where / = 1 is assumed. In our case there 
exist up to 3 CPs, and hence we have to extend the Lam- 
bert W function approach as follows. First, when 3 CPs 
exist, only the innermost and outermost ones turn out 
relevant. We evaluate D by choosing each of them, one 
after another, as the SP. In some portion of the com- 
putational domain (y > 1), D for one CP may exceed 
1/e, and hence the solution is not defined. Call this solu- 
tion the "broken solution" , denoted by Mb. Choosing the 
other CP as SP results in a continuous solution, denoted 
by Mc- If standing shocks exist, they have to appear 
where the Rankine-Hugoniot relations and evolutionary 
condition are met. In the isothermal case these translate 
into [4, 12] 



M+M- = 1 and M+ > 1, 



(5) 



respectively. Here -I- (— ) represents the shock upstream 
(downstream). This suggests a simple graphical means 
to construct solutions with shocks [l2[, where we plot Mb, 
and examine whether it intersects the 1/Mc curve. Any 
intersection represents a shock jump, however the solu- 
tion cannot jump from a lower to a higher curve. 

Figure [2] illustrates our solution procedure, giving the 
radial dependence of the Mach number M ((a) and (c)) 
and D ((b) and (d)). In FigsUb and Hi, the light hor- 
izontal lines represent 1/e. The solid and dashed lines 




FIG. 2: Wind solutions with the streamer geometry for which 
/oo = 6, /m = 14, rc = 3i?o, and 5 = 0.7 Rq. In the left 
and right columns, T is 1.2 and 1.3 MK, respectively. Panels 
(a) and (c) present the r— dependence of the Mach number 
M, while (b) and (d) give that of D. The solid (dashed) lines 
represent the continuous (broken) solutions. In (b) and (d), 
the horizontal lines give 1 /e for comparison. Standing shocks 
are allowed only when the two curves 1/Mc and Mi, intersect. 



correspond to the continuous and broken solutions, re- 
spectively. In addition to M, Figsl^^ and also give 
1/Mc. Figures Hi and [2}d are for T = 1.2 MK, while 
FigsHt and Hi are for T = 1.3 MK. In both cases the 
tube parameters are foo = 6, /a/ = 14, rc = 3i?0, and 
6 = 0.7 Rq. Consider now FigsH^ andH". It is seen that 
both curves in Figl2}D exhibit three local extrema, whose 
locations correspond to the CPs. This follows from that 
dD/dy = at any CP (see Eq.dH)). Furthermore, the 
global maximum of D is attained at the outermost CP, 
located at 4.89i?o. Therefore when the innermost CP is 
chosen as the SP, D > 1/e around the outermost CP for 
4.2 < r < 6.51i?o. RecaUing that W{-D) is real-valued 
only when —D > —1/e, one readily understands that in 
this interval choosing the innermost CP as the SP does 
not result in a solution to Eq.(IT]). Figure [5^ also shows 
that the curve 1/Mc does not intersect Mb, indicating the 
solution to Eq.([T]) is unique and is the continuous one. 

The situation changes when T = 1.3 MK. Now the 
global maximum of D is attained at the innermost CP, 
located at 1.75i?o (Fig Hi). Choosing the outermost CP 
as the SP leads to that D > 1/e in the interval [1.53, 
1.98]i?o where there is no solution (FigH;). However, 
two standing shocks are now allowed, since two crossings 
exist between the curves 1/Mc and Mb, located at 2.11 
and 3.96i?o, respectively. Hence in addition to the con- 
tinuous one (Mc adopting the innermost CP as the SP), 
two additional solutions exist to Eq.([T|): both start with 
Mc but one connects to Mb at the inner crossing, the 
other connects to Mb at the outer one. 

Although Eq.(IT|) permits solutions with shocks, and 
time-dependent simulations suggest these steady-state 
solutions can be attained[l,|3,[l,[3], one may still question 
whether the shock solutions can stand the sensitivity test 
similar to@ which showed standing shocks in the solar 
wind from the center of coronal holes are very unlikely. 
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FIG. 3: Regions in the T — rc space where standing shocks 
are allowed, shown as the area bounded by the two curves in 
same color. The results are obtained by varying fao, /m and 
5 respectively, the rest of the parameters are fixed at their 
reference values /oo = 6, /m = 14, 5 = 0.7 Rq. 



for the parameter range allowing shock solutions is ex- 
tremely limited. To see whether the same happens with 
the streamer geometry, we note that given /oo, Jm and 
(5, shock solutions are allowed only in the area bounded 
by two curves in the [T, rc] space. Figure [3] presents a 
series of such curves obtained by varying (a) /oo, (b) /a/, 
and (c) S about the reference values /oo = 6, /m — 14, 
and S = O.7i?0. (In what follows, the temperatures are 
in MK, rc and d in Rq.) Let us first examine the cases 
with reference values (the black curves connecting as- 
terisks) . Figure [3] shows that the area bounded by the 
two curves is rather broad, and with increasing rc, both 
curves are shifted towards lower temperatures, indicating 
streamers whose cusps are located higher in the corona 
are more likely associated with standing shocks. For in- 
stance, when the cusp is located at 3.6i?0, the slow wind 
may possess standing shocks as long as 0.91 <T < 2.04, 
which actually tends to be lower than the often-quoted 
values of coronal temperatures. On the other hand, even 
for the lowest cusp height examined {rc = 2.4Rq), the 
temperature range is [1.74,2.80], still largely compatible 
with the observational range. 

Figure [3^ examines the effects of varying /oo- It is 
seen that while /oo decreases from its reference value to 
2 (the red curves), the T range allowing shocks increases 
significantly. Actually for rc > 2.5, in the examined tem- 
perature range there virtually exists no upper bound for 
shocks to occur. Take rc = 3.6 for instance. Shock so- 
lutions take place as long as T > 1.35. On the contrary, 
increasing /oo to 10 (the blue curves) makes shocks ap- 
pear in a much narrower T range (the width is 0.3). 



The effects of varying /m are shown in FiglSlD, which 
shows that increasing /m considerably broadens the area 
allowing standing shocks. For example, with /m increas- 
ing from 6 to 14, the width along the T-axis of the area 
increases from ~ 0.1 to 1.1. When /m further increases 
to 22, this width increases dramatically from 0.88 at 
rc = 2.4 to > 2.9 for r^ > 2.6. Figure Efc shows what 
happens when 6 changes, where it is seen that increasing 
6 reduces the range of T where shocks are allowed. For 
instance, with rc = 3.0, this range for S = 0.4 {6 = 1) 
is [1.27,3.21] ([1.55,2.26]), while the range for the ref- 
erence value S = 0.7 lies in between. It is interesting 
to note that for S = I, at rc ^ 2.68 the upper bound 
for T (the right blue curve) changes its slope dramati- 
cally, and for rc < 2.53 no shock solutions exist. For 
2.53 < rc < 2.68, it turns out that on the right of the 
right blue curve actually no solution exists, since now 
only two critical points exist and neither of them cor- 
responds to a Z? < 1/e throughout the computational 
domain (see Eq.®). This is different from the portion 
rc > 2.68, where on the right of the right blue curve 
there does exist a solution which is the continuous one. 
Putting the three panels together, one may see that for 
most combinations of tube parameters, the area in the 
T ~ rc space supporting standing shocks is substantial. 
Hence with the streamer geometry, standing shocks in 
the inner slow wind seem physically accessible. 

It is not easy to exhaust the possible tube parameters 
and the consequent changes in shock properties. Let us 
instead discuss only the shocks found, examining their 
detectability. First, Sp, the density jump relative to the 
upstream value, is up to 8, a result of the isothermal 
assumption exceeding the nominal upper limit of 4 for 
adiabatic gases. As shown by[l3], a. Sp oi ^ 2.3 at a 
standing shock produces an enhancement in the polarized 
brightness intensity that is only marginally detectable. A 
dp of 8 certainly makes such detections easier, but one 
can not say this for sure without constructing detailed 
observables. Second, by conserving angular momentum 
a coronal shock also produces a discontinuity in the az- 
imuthal fiow speed v^, leading in principle to measurable 
Doppler shifts in H I Ly a. However, the jump in 
turns out < 4 km/s, discerning which is way beyond the 
sensitivity of SOHO/UVCS, whose spectral resolution of 
0.23 A translates into ^ 57 km/s. 

The isothermal assumption needs some justification. 
First, it is not far from reality. The UVCS measurements 
of the H I Ly a emission from an equatorial streamer [TJ] 
showed that the proton kinetic temperature Tp in the 
stalk decreases only mildly from 1.45 MK at 3.6i?0 to 
1.3 MK at 5.1i?0 (their Fig.3b). If the stalk and one 
of streamer legs are on the same flow tube, then Fig. 4b 
injlll shows that 1.41 < Tp < 2.09 MK at 2.33i?0 (the 
leftmost two open circles and rightmost two solid ones 
in their Fig.4d). As for Tg, the electron-scattered H I 
Lya measured by UVCS yielded a of 1.1 ± 0.3 MK 
at 2.7i?0[lH|- Although for a streamer, this value may 
serve to estimate Tg in flowing regions at similar heights. 
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Direct Tg measurements above that distance are sparse. 
Nonetheless, multi-fluid MHD models indicate that Tg 
ranges from 0.8 MK at 3i?0 to 0.65 MK at 5Rq (Fig.Sd 
in [17]). The mean of and Tp, the temperature T in 
this study is thus ^ 1.1 — 1.8 MK at 2.3Rq and decreases 
to ~ IMK at 5i?0. Furthermore, T at the slow wind 
source region is ^ 0.8 — 1.2 MK, be this source in a coro- 
nal hole or in its neighboring quiet Sun [16]. Second, in- 



troducing a more complete energy equation, as was done 
in|5| for a coronal-hole how, will likely strengthen rather 
than weaken our conclusion. That study shows that in- 
troducing thermal conduction and two- fluid effects allows 
for a much broader parameter range supporting standing 
shocks, compared with isothermal and polytropic com- 
putations. 
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